Adaptive optical third-harmonic generation microscopy for in vivo imaging of tissues.

Third-harmonic generation microscopy is a powerful label-free nonlinear imaging technique, providing essential information about structural characteristics of cells and tissues without requiring external labelling agents. In this work, we integrated a recently developed compact adaptive optics module into a third-harmonic generation microscope, to measure and correct for optical aberrations in complex tissues. Taking advantage of the high sensitivity of the third-harmonic generation process to material interfaces and thin membranes, along with the 1,300-nm excitation wavelength used here, our adaptive optical third-harmonic generation microscope enabled high-resolution in vivo imaging within highly scattering biological model systems. Examples include imaging of myelinated axons and vascular structures within the mouse spinal cord and deep cortical layers of the mouse brain, along with imaging of key anatomical features in the roots of the model plant Brachypodium distachyon . In all instances, aberration correction led to significant enhancements in image quality.

CO2 Response Screen in Grass Brachypodium Reveals Key Role of a MAP Kinase in CO2-Triggered Stomatal Closure.

Plants respond to increased CO2 concentrations through stomatal closure, which can contribute to increased water use efficiency. Grasses display faster stomatal responses than eudicots due to dumbbell-shaped guard cells flanked by subsidiary cells working in opposition. However, forward genetic screening for stomatal CO2 signal transduction mutants in grasses has yet to be reported. The grass model Brachypodium distachyon is closely related to agronomically important cereal crops, sharing largely collinear genomes. To gain insights into CO2 control mechanisms of stomatal movements in grasses, we developed an unbiased forward genetic screen with an EMS-mutagenized Brachypodium distachyon M5 generation population using infrared imaging to identify plants with altered leaf temperatures at elevated CO2. Among isolated mutants, a "chill1" mutant exhibited cooler leaf temperatures than wildtype Bd21-3 parent control plants after exposure to increased [CO2]. chill1 plants showed strongly impaired high CO2-induced stomatal closure despite retaining a robust abscisic acid-induced stomatal closing response. Through bulked segregant whole-genome-sequencing analyses followed by analyses of further backcrossed F4 generation plants and generation and characterization of sodium-azide and CRISPR-cas9 mutants, chill1 was mapped to a protein kinase, Mitogen-Activated Protein Kinase 5 (BdMPK5). The chill1 mutation impaired BdMPK5 protein-mediated CO2/HCO3- sensing together with the High Temperature 1 (HT1) Raf-like kinase in vitro. Furthermore, AlphaFold2-directed structural modeling predicted that the identified BdMPK5-D90N chill1 mutant residue is located at the interface of BdMPK5 with the BdHT1 Raf-like kinase. BdMPK5 is a key signaling component that mediates CO2-induced stomatal movements and is proposed to function as a component of the primary CO2 sensor in grasses.

Transposition of HOPPLA in siRNA-deficient plants suggests a limited effect of the environment on retrotransposon mobility in Brachypodium distachyon.

Long terminal repeat retrotransposons (LTR-RTs) are powerful mutagens regarded as a major source of genetic novelty and important drivers of evolution. Yet, the uncontrolled and potentially selfish proliferation of LTR-RTs can lead to deleterious mutations and genome instability, with large fitness costs for their host. While population genomics data suggest that an ongoing LTR-RT mobility is common in many species, the understanding of their dual role in evolution is limited. Here, we harness the genetic diversity of 320 sequenced natural accessions of the Mediterranean grass Brachypodium distachyon to characterize how genetic and environmental factors influence plant LTR-RT dynamics in the wild. When combining a coverage-based approach to estimate global LTR-RT copy number variations with mobilome-sequencing of nine accessions exposed to eight different stresses, we find little evidence for a major role of environmental factors in LTR-RT accumulations in B. distachyon natural accessions. Instead, we show that loss of RNA polymerase IV (Pol IV), which mediates RNA-directed DNA methylation in plants, results in high transcriptional and transpositional activities of RLC_BdisC024 (HOPPLA) LTR-RT family elements, and that these effects are not stress-specific. This work supports findings indicating an ongoing mobility in B. distachyon and reveals that host RNA-directed DNA methylation rather than environmental factors controls their mobility in this wild grass model.

The reference genome and abiotic stress responses of the model perennial grass Brachypodium sylvaticum.

Perennial grasses are important forage crops and emerging biomass crops and have the potential to be more sustainable grain crops. However, most perennial grass crops are difficult experimental subjects due to their large size, difficult genetics, and/or their recalcitrance to transformation. Thus, a tractable model perennial grass could be used to rapidly make discoveries that can be translated to perennial grass crops. Brachypodium sylvaticum has the potential to serve as such a model because of its small size, rapid generation time, simple genetics, and transformability. Here, we provide a high-quality genome assembly and annotation for B. sylvaticum, an essential resource for a modern model system. In addition, we conducted transcriptomic studies under 4 abiotic stresses (water, heat, salt, and freezing). Our results indicate that crowns are more responsive to freezing than leaves which may help them overwinter. We observed extensive transcriptional responses with varying temporal dynamics to all abiotic stresses, including classic heat-responsive genes. These results can be used to form testable hypotheses about how perennial grasses respond to these stresses. Taken together, these results will allow B. sylvaticum to serve as a truly tractable perennial model system.

JGI Plant Gene Atlas: an updateable transcriptome resource to improve functional gene descriptions across the plant kingdom.

Gene functional descriptions offer a crucial line of evidence for candidate genes underlying trait variation. Conversely, plant responses to environmental cues represent important resources to decipher gene function and subsequently provide molecular targets for plant improvement through gene editing. However, biological roles of large proportions of genes across the plant phylogeny are poorly annotated. Here we describe the Joint Genome Institute (JGI) Plant Gene Atlas, an updateable data resource consisting of transcript abundance assays spanning 18 diverse species. To integrate across these diverse genotypes, we analyzed expression profiles, built gene clusters that exhibited tissue/condition specific expression, and tested for transcriptional response to environmental queues. We discovered extensive phylogenetically constrained and condition-specific expression profiles for genes without any previously documented functional annotation. Such conserved expression patterns and tightly co-expressed gene clusters let us assign expression derived additional biological information to 64 495 genes with otherwise unknown functions. The ever-expanding Gene Atlas resource is available at JGI Plant Gene Atlas (https://plantgeneatlas.jgi.doe.gov) and Phytozome (https://phytozome.jgi.doe.gov/), providing bulk access to data and user-specified queries of gene sets. Combined, these web interfaces let users access differentially expressed genes, track orthologs across the Gene Atlas plants, graphically represent co-expressed genes, and visualize gene ontology and pathway enrichments.

PHYTOCHROME C regulation of photoperiodic flowering via PHOTOPERIOD1 is mediated by EARLY FLOWERING 3 in Brachypodium distachyon.

Daylength sensing in many plants is critical for coordinating the timing of flowering with the appropriate season. Temperate climate-adapted grasses such as Brachypodium distachyon flower during the spring when days are becoming longer. The photoreceptor PHYTOCHROME C is essential for long-day (LD) flowering in B. distachyon. PHYC is required for the LD activation of a suite of genes in the photoperiod pathway including PHOTOPERIOD1 (PPD1) that, in turn, result in the activation of FLOWERING LOCUS T (FT1)/FLORIGEN, which causes flowering. Thus, B. distachyon phyC mutants are extremely delayed in flowering. Here we show that PHYC-mediated activation of PPD1 occurs via EARLY FLOWERING 3 (ELF3), a component of the evening complex in the circadian clock. The extreme delay of flowering of the phyC mutant disappears when combined with an elf3 loss-of-function mutation. Moreover, the dampened PPD1 expression in phyC mutant plants is elevated in phyC/elf3 mutant plants consistent with the rapid flowering of the double mutant. We show that loss of PPD1 function also results in reduced FT1 expression and extremely delayed flowering consistent with results from wheat and barley. Additionally, elf3 mutant plants have elevated expression levels of PPD1, and we show that overexpression of ELF3 results in delayed flowering associated with a reduction of PPD1 and FT1 expression, indicating that ELF3 represses PPD1 transcription consistent with previous studies showing that ELF3 binds to the PPD1 promoter. Indeed, PPD1 is the main target of ELF3-mediated flowering as elf3/ppd1 double mutant plants are delayed flowering. Our results indicate that ELF3 operates downstream from PHYC and acts as a repressor of PPD1 in the photoperiod flowering pathway of B. distachyon.

Mixed-Linkage Glucan Is the Main Carbohydrate Source and Starch Is an Alternative Source during Brachypodium Grain Germination.

Seeds of the model grass Brachypodium distachyon are unusual because they contain very little starch and high levels of mixed-linkage glucan (MLG) accumulated in thick cell walls. It was suggested that MLG might supplement starch as a storage carbohydrate and may be mobilised during germination. In this work, we observed massive degradation of MLG during germination in both endosperm and nucellar epidermis. The enzymes responsible for the MLG degradation were identified in germinated grains and characterized using heterologous expression. By using mutants targeting MLG biosynthesis genes, we showed that the expression level of genes coding for MLG and starch-degrading enzymes was modified in the germinated grains of knocked-out cslf6 mutants depleted in MLG but with higher starch content. Our results suggest a substrate-dependent regulation of the storage sugars during germination. These overall results demonstrated the function of MLG as the main carbohydrate source during germination of Brachypodium grain. More astonishingly, cslf6 Brachypodium mutants are able to adapt their metabolism to the lack of MLG by modifying the energy source for germination and the expression of genes dedicated for its use.

Regulation of hair cell and stomatal size by a hair cell-specific peroxidase in the grass Brachypodium distachyon.

The leaf epidermis is the outermost cell layer forming the interface between plants and the atmosphere that must both provide a robust barrier against (a)biotic stressors and facilitate carbon dioxide uptake and leaf transpiration.1 To achieve these opposing requirements, the plant epidermis developed a wide range of specialized cell types such as stomata and hair cells. Although factors forming these individual cell types are known,2,3,4,5 it is poorly understood how their number and size are coordinated. Here, we identified a role for BdPRX76/BdPOX, a class III peroxidase, in regulating hair cell and stomatal size in the model grass Brachypodium distachyon. In bdpox mutants, prickle hair cells were smaller and stomata were longer. Because stomatal density remained unchanged, the negative correlation between stomatal size and density was disrupted in bdpox and resulted in higher stomatal conductance and lower intrinsic water-use efficiency. BdPOX was exclusively expressed in hair cells, suggesting that BdPOX cell-autonomously promotes hair cell size and indirectly restricts stomatal length. Cell-wall autofluorescence and lignin stainings indicated a role for BdPOX in the lignification or crosslinking of related phenolic compounds at the hair cell base. Ectopic expression of BdPOX in the stomatal lineage increased phenolic autofluorescence in guard cell (GC) walls and restricted stomatal elongation in bdpox. Together, we highlight a developmental interplay between hair cells and stomata that optimizes epidermal functionality. We propose that cell-type-specific changes disrupt this interplay and lead to compensatory developmental defects in other epidermal cell types.

Multiple origins, one evolutionary trajectory: gradual evolution characterizes distinct lineages of allotetraploid Brachypodium.

The "genomic shock" hypothesis posits that unusual challenges to genome integrity such as whole genome duplication may induce chaotic genome restructuring. Decades of research on polyploid genomes have revealed that this is often, but not always the case. While some polyploids show major chromosomal rearrangements and derepression of transposable elements in the immediate aftermath of whole genome duplication, others do not. Nonetheless, all polyploids show gradual diploidization over evolutionary time. To evaluate these hypotheses, we produced a chromosome-scale reference genome for the natural allotetraploid grass Brachypodium hybridum, accession "Bhyb26." We compared 2 independently derived accessions of B. hybridum and their deeply diverged diploid progenitor species Brachypodium stacei and Brachypodium distachyon. The 2 B. hybridum lineages provide a natural timecourse in genome evolution because one formed 1.4 million years ago, and the other formed 140 thousand years ago. The genome of the older lineage reveals signs of gradual post-whole genome duplication genome evolution including minor gene loss and genome rearrangement that are missing from the younger lineage. In neither B. hybridum lineage do we find signs of homeologous recombination or pronounced transposable element activation, though we find evidence supporting steady post-whole genome duplication transposable element activity in the older lineage. Gene loss in the older lineage was slightly biased toward 1 subgenome, but genome dominance was not observed at the transcriptomic level. We propose that relaxed selection, rather than an abrupt genomic shock, drives evolutionary novelty in B. hybridum, and that the progenitor species' similarity in transposable element load may account for the subtlety of the observed genome dominance.

Brachypodium: 20 years as a grass biology model system; the way forward?

It has been 20 years since Brachypodium distachyon was suggested as a model grass species, but ongoing research now encompasses the entire genus. Extensive Brachypodium genome sequencing programmes have provided resources to explore the determinants and drivers of population diversity. This has been accompanied by cytomolecular studies to make Brachypodium a platform to investigate speciation, polyploidisation, perenniality, and various aspects of chromosome and interphase nucleus organisation. The value of Brachypodium as a functional genomic platform has been underscored by the identification of key genes for development, biotic and abiotic stress, and cell wall structure and function. While Brachypodium is relevant to the biofuel industry, its impact goes far beyond that as an intriguing model to study climate change and combinatorial stress.

N-dependent dynamics of root growth and nitrate and ammonium uptake are altered by the bacterium Herbaspirillum seropedicae in the cereal model Brachypodium distachyon.

Nitrogen (N) fixation in cereals by root-associated bacteria is a promising solution for reducing use of chemical N fertilizers in agriculture. However, plant and bacterial responses are unpredictable across environments. We hypothesized that cereal responses to N-fixing bacteria are dynamic, depending on N supply and time. To quantify the dynamics, a gnotobiotic, fabricated ecosystem (EcoFAB) was adapted to analyse N mass balance, to image shoot and root growth, and to measure gene expression of Brachypodium distachyon inoculated with the N-fixing bacterium Herbaspirillum seropedicae. Phenotyping throughput of EcoFAB-N was 25-30 plants h-1 with open software and imaging systems. Herbaspirillum seropedicae inoculation of B. distachyon shifted root and shoot growth, nitrate versus ammonium uptake, and gene expression with time; directions and magnitude depended on N availability. Primary roots were longer and root hairs shorter regardless of N, with stronger changes at low N. At higher N, H. seropedicae provided 11% of the total plant N that came from sources other than the seed or the nutrient solution. The time-resolved phenotypic and molecular data point to distinct modes of action: at 5 mM NH4NO3 the benefit appears through N fixation, while at 0.5 mM NH4NO3 the mechanism appears to be plant physiological, with H. seropedicae promoting uptake of N from the root medium.Future work could fine-tune plant and root-associated microorganisms to growth and nutrient dynamics.

The recent evolutionary rescue of a staple crop depended on over half a century of global germplasm exchange.

Rapid environmental change can lead to population extinction or evolutionary rescue. The global staple crop sorghum (Sorghum bicolor) has recently been threatened by a global outbreak of an aggressive new biotype of sugarcane aphid (SCA; Melanaphis sacchari). We characterized genomic signatures of adaptation in a Haitian breeding population that had rapidly adapted to SCA infestation, conducting evolutionary population genomics analyses on 296 Haitian lines versus 767 global accessions. Genome scans and geographic analyses suggest that SCA adaptation has been conferred by a globally rare East African allele of RMES1, which spread to breeding programs in Africa, Asia, and the Americas. De novo genome sequencing revealed potential causative variants at RMES1. Markers developed from the RMES1 sweep predicted resistance in eight independent commercial and public breeding programs. These findings demonstrate the value of evolutionary genomics to develop adaptive trait technology and highlight the benefits of global germplasm exchange to facilitate evolutionary rescue.

Cell Wall Compositions of Sorghum bicolor Leaves and Roots Remain Relatively Constant Under Drought Conditions.

Renewable fuels are needed to replace fossil fuels in the immediate future. Lignocellulosic bioenergy crops provide a renewable alternative that sequesters atmospheric carbon. To prevent displacement of food crops, it would be advantageous to grow biofuel crops on marginal lands. These lands will likely face more frequent and extreme drought conditions than conventional agricultural land, so it is crucial to see how proposed bioenergy crops fare under these conditions and how that may affect lignocellulosic biomass composition and saccharification properties. We found that while drought impacts the plant cell wall of Sorghum bicolor differently according to tissue and timing of drought induction, drought-induced cell wall compositional modifications are relatively minor and produce no negative effect on biomass conversion. This contrasts with the cell wall-related transcriptome, which had a varied range of highly variable genes (HVGs) within four cell wall-related GO categories, depending on the tissues surveyed and time of drought induction. Further, many HVGs had expression changes in which putative impacts were not seen in the physical cell wall or which were in opposition to their putative impacts. Interestingly, most pre-flowering drought-induced cell wall changes occurred in the leaf, with matrix and lignin compositional changes that did not persist after recovery from drought. Most measurable physical post-flowering cell wall changes occurred in the root, affecting mainly polysaccharide composition and cross-linking. This study couples transcriptomics to cell wall chemical analyses of a C4 grass experiencing progressive and differing drought stresses in the field. As such, we can analyze the cell wall-specific response to agriculturally relevant drought stresses on the transcriptomic level and see whether those changes translate to compositional or biomass conversion differences. Our results bolster the conclusion that drought stress does not substantially affect the cell wall composition of specific aerial and subterranean biomass nor impede enzymatic hydrolysis of leaf biomass, a positive result for biorefinery processes. Coupled with previously reported results on the root microbiome and rhizosphere and whole transcriptome analyses of this study, we can formulate and test hypotheses on individual gene candidates' function in mediating drought stress in the grass cell wall, as demonstrated in sorghum.

Plant single-cell solutions for energy and the environment.

Progress in sequencing, microfluidics, and analysis strategies has revolutionized the granularity at which multicellular organisms can be studied. In particular, single-cell transcriptomics has led to fundamental new insights into animal biology, such as the discovery of new cell types and cell type-specific disease processes. However, the application of single-cell approaches to plants, fungi, algae, or bacteria (environmental organisms) has been far more limited, largely due to the challenges posed by polysaccharide walls surrounding these species' cells. In this perspective, we discuss opportunities afforded by single-cell technologies for energy and environmental science and grand challenges that must be tackled to apply these approaches to plants, fungi and algae. We highlight the need to develop better and more comprehensive single-cell technologies, analysis and visualization tools, and tissue preparation methods. We advocate for the creation of a centralized, open-access database to house plant single-cell data. Finally, we consider how such efforts should balance the need for deep characterization of select model species while still capturing the diversity in the plant kingdom. Investments into the development of methods, their application to relevant species, and the creation of resources to support data dissemination will enable groundbreaking insights to propel energy and environmental science forward.

Identification of Candidate Susceptibility Genes to Puccinia graminis f. sp. tritici in Wheat.

Wheat stem rust disease caused by Puccinia graminis f. sp. tritici (Pgt) is a global threat to wheat production. Fast evolving populations of Pgt limit the efficacy of plant genetic resistance and constrain disease management strategies. Understanding molecular mechanisms that lead to rust infection and disease susceptibility could deliver novel strategies to deploy crop resistance through genetic loss of disease susceptibility. We used comparative transcriptome-based and orthology-guided approaches to characterize gene expression changes associated with Pgt infection in susceptible and resistant Triticum aestivum genotypes as well as the non-host Brachypodium distachyon. We targeted our analysis to genes with differential expression in T. aestivum and genes suppressed or not affected in B. distachyon and report several processes potentially linked to susceptibility to Pgt, such as cell death suppression and impairment of photosynthesis. We complemented our approach with a gene co-expression network analysis to identify wheat targets to deliver resistance to Pgt through removal or modification of putative susceptibility genes.

Plant Pan-Genomics Comes of Age.

A pan-genome is the nonredundant collection of genes and/or DNA sequences in a species. Numerous studies have shown that plant pan-genomes are typically much larger than the genome of any individual and that a sizable fraction of the genes in any individual are present in only some genomes. The construction and interpretation of plant pan-genomes are challenging due to the large size and repetitive content of plant genomes. Most pan-genomes are largely focused on nontransposable element protein coding genes because they are more easily analyzed and defined than noncoding and repetitive sequences. Nevertheless, noncoding and repetitive DNA play important roles in determining the phenotype and genome evolution. Fortunately, it is now feasible to make multiple high-quality genomes that can be used to construct high-resolution pan-genomes that capture all the variation. However, assembling, displaying, and interacting with such high-resolution pan-genomes will require the development of new tools.

Gradual polyploid genome evolution revealed by pan-genomic analysis of Brachypodium hybridum and its diploid progenitors.

Our understanding of polyploid genome evolution is constrained because we cannot know the exact founders of a particular polyploid. To differentiate between founder effects and post polyploidization evolution, we use a pan-genomic approach to study the allotetraploid Brachypodium hybridum and its diploid progenitors. Comparative analysis suggests that most B. hybridum whole gene presence/absence variation is part of the standing variation in its diploid progenitors. Analysis of nuclear single nucleotide variants, plastomes and k-mers associated with retrotransposons reveals two independent origins for B. hybridum, ~1.4 and ~0.14 million years ago. Examination of gene expression in the younger B. hybridum lineage reveals no bias in overall subgenome expression. Our results are consistent with a gradual accumulation of genomic changes after polyploidization and a lack of subgenome expression dominance. Significantly, if we did not use a pan-genomic approach, we would grossly overestimate the number of genomic changes attributable to post polyploidization evolution.

The fate of 35S rRNA genes in the allotetraploid grass Brachypodium hybridum.

Nucleolar dominance (ND) consists of the reversible silencing of 35S/45S rDNA loci inherited from one of the ancestors of an allopolyploid. The molecular mechanisms by which one ancestral rDNA set is selected for silencing remain unclear. We applied a combination of molecular (Southern blot hybridization and reverse-transcription cleaved amplified polymorphic sequence analysis), genomic (analysis of variants) and cytogenetic (fluorescence in situ hybridization) approaches to study the structure, expression and epigenetic landscape of 35S rDNA in an allotetraploid grass that exhibits ND, Brachypodium hybridum (genome composition DDSS), and its putative progenitors, Brachypodium distachyon (DD) and Brachypodium stacei (SS). In progenitor genomes, B. stacei showed a higher intragenomic heterogeneity of rDNA compared with B. distachyon. In all studied accessions of B. hybridum, there was a reduction in the copy number of S homoeologues, which was accompanied by their inactive transcriptional status. The involvement of DNA methylation in CG and CHG contexts in the silencing of the S-genome rDNA loci was revealed. In the B. hybridum allotetraploid, ND is stabilized towards the D-genome units, irrespective of the polyphyletic origin of the species, and does not seem to be influenced by homoeologous 35S rDNA ratios and developmental stage.